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We suggest a procedure for demonstrating quantum coherence and measuring decoherence times 
between different fluxoid states of a SQUID by using "adiabatic inversion", where one macroscopic 
fluxoid state is smoothly transferred into the other, like a spin reversing direction by following 
a slowly moving magnetic field. This is accomplished by sweeping an external applied flux, and 
depends on a well-deflned quantum phase between the two macroscopic states. Varying the speed of 
the sweep relative to the decoherence time permits one to move from the quantum regime, where such 
a well-defined phase exists, to the classical regime where it is lost and the inversion is inhibited. Thus 
observing whether inversion has taken place or not as a function of sweep speed offers the possibility 
of measuring the decoherence time. The main requirement for the feasibility of the scheme appears 
to be that the low temperature relaxation time among the quantum levels of the SQUID be long 
compared to other time scales of the problem. Applications to the "quantum computer", with the 
level system of the SQUID playing the role of the qbit, are briefiy discussed. 



I. INTRODUCTION 

One of the most interesting questions of "macro- 
scopic quantum mechanics" concerns the demonstration 
of quantum coherence and interference effects between 
apparently grossly distinguishable states of large systems. 
Observation of such effects would definitively lay to rest 
any ideas about big or "classical" systems being funda- 
mentally different than small systems, in some way not 
subject to the rules of quantum mechanics. One system 
which has been extensively discussed in this regard is the 
Josephson junction ^ under current bias and the closely 
related rf Squid, where a superconducting ring is inter- 
rupted by a Josephson junction. In these systems the 
phase across the junction or the flux $ in the ring plays 
the role of a collective coordinate. In the last decade a 
number of beautiful junction experiments at low temper- 
ature Q have seen effects connected with the quantized 
energy levels Q expected with respect to this coordinate. 
A fast sweeping method [Q has also seen the effects of 
these quantized levels even at relatively high tempera- 
ture. 

An even more striking and direct manifestation of 
macroscopic quantum behavior would be the demon- 
stration of quantum coherence between two apparently 
macroscopically different states of the rf SQUID, between 
states where the current goes around the ring in oppo- 
site directions. Recently, results have been reported with 
microwave methods showing the repulsion of energy lev- 
els expected due to quantum mechanical mixing of the 



different fluxoid states [^|. Here we would like to sug- 
gest a new method for demonstrating macroscopic co- 
herence, one which further allows a measurement of the 
decoherence time, the time in which the coherence be- 
tween the states is lost. Our basic idea is to use the pro- 
cess of "adiabatic inversion" where a slowly varying held 
is used to reverse the states of a quantum system. In its 
most straightforward realization, our proposal consists of 
starting with the system in its lowest state, making a fast 
but adiabatic sweep, and reading out to see if the flnal 
state is the same or the opposite fluxoid state. If the state 
has switched, the system has behaved quantum mechan- 
ically, with phase coherence between the two states. If 
not, the phase coherence between the states was lost and 
the system stayed in its original conflguration- behaved 
classically. 

A very interesting aspect of the present proposal, as 
we shall discuss below, is the possiblity of passing from 
the quantum to the classical regime by simply varying 
the sweep speed. This allows us to obtain, for adiabatic 
conditions, the decoherence time as the longest sweep 
time for which the inversion is successful. 



II. ADIABATIC INVERSION IN THE RF SQUID 

For the rf SQUID biased with an external flux $2,, 
the "potential energy" of the system for which $ is the 
coordinate is shown in Fig 1. One has, near the origin, 
the familiar double well potential 
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C/ = C/o[l/2($-$,)2 + /3icos$] . (1) 

Depending on whether the system is in the left or right 
well, the flux through the ring has a different sign and the 
current goes around in opposite directions. The poten- 
tial can be varied by altering the external ^x, becoming 
symmetric for $j; = 0. We suppose a situation where 
a few states may be possible in a well, but we will be 
essentially concerned with the lowest level in each. In 
the adiabatic inversion procedure to be discussed, is 
swept from a maximum to an minimum value, passing 
through zero, such that the initially asymmetric left and 
right wells exchange roles, the originally higher well be- 
coming the lower one and vice versa. The asymmetry 
of the configurations, however, is kept small so that we 
effectively have only a two-state system, composed of the 
lowest state in each well. 




FIG. 1. Double potential well with harmonic level spacing 
uio and initial spacing e between lowest levels. 

As quantum states we refer to the lowest state of the 
left or right well as |i > and \R >, and the tunneling 
between \L > and \R > will generally lead to eigenstates 
which are linear combinations of the two. Any two-state 
system may viewed as constituting the two components 
of a "spin" , and this provides an easy visualization which 
has been used in many contexts @] . In the present prob- 
lem "spin up" can be identified with, say, the state \L > 
(flux one way) and "spin down" with \R > (flux the other 
way), while the spin in a general direction represents a 
quantum mechanical linear combination of the two states, 
with some definite relative magnitude and phase. Various 
influences like the external bias or the tunneling ampli- 
tude may be thought of as creating a kind of pseudo- 
magnetic field V about which the spin or polarization P 
precesses (Fig. 2) according to the equation 

P = V X P - DPtr , (2) 

where V can be time dependent. The quantity D is the 
decoherence parameter, which we neglect for the moment 
and will deal with below. 

Now it is a familiar fact under adiabatic conditions, 
where V varies slowly, that the "spin" P will tend to 
"follow" the magnetic field V(t). This is a completely 



familiar procedure when rotating the spin of say an atom 
or a neutron by a magnetic field. Here the motion of V 
will be accomplished by sweeping ^x', it is then merely 
necessary to identify the components of V and the mean- 
ing of "adiabatic" or "slow" . 




FIG. 2. Precession of the "spin" vector P around the 
pseudo-magnetic field vector V(t). In adiabatic inversion V, 
swings from "up" to "down", and carries P with it. 



The vertical component of V corresponds to the dif- 
ference in the two potential wells or lowest energy lev- 
els, for small asymmetry it is linear in ^x- If e is the 
initial level splitting when ^x — we can write 

Vvertit) = e($a;(t)/<l)™''^). Thus as ^x sweeps from its 
positive maximum value to its negative minimum value 
Vuert reverses direction. The transverse component of 
'Vtr corresponds to the tunneling energy between the two 
quasi-degenerate states Vtr = ^tunnel- As Vyert passes 
through zero aX ^x — the \L > and \R > states are 
strongly mixed and the splitting of the resulting energy 
eigenstates is determined by Vtr alone. The magnitude 
of V{t) at a given time, |V| = y/V^^^t + Vtl 

gives the in- 
stantaneous splitting of the two levels, which varies from 
approximately e in the vertical position of V to LOtunnei 
in the horizontal position. 

As to the meaning of "adiabatic" or "slow", the im- 
posed time variation of V should not significantly con- 
tain frequencies or fourier components corresponding to 
the energy splitting between levels. Or in terms of time, 
the rate of variation of V should be slow on the time 
scale corresponding to the tunneling time between the 
two states. Thus we have the requirement on V that 
its relative rate of variation V /V always be small com- 
pared to V itself. Since the varying component of V is 
Vverticai (neglecting the indirect effect of ^x on ujtunnei) 
we require Vverticai /V « V. Thus taking the near de- 
generate configuration where V w totunnei we find 

^x{t) _ 2 /ON 
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as the condition for adiabaticity. If the adiabatic condi- 
tion is violated, then P cannot follow V and undergoes 
wide swings. We stress that for adiabatic inversion there 
must be a well-defined quantum phase between the two 
states, and that when this phase is lost the inversion is 
suppressed. This is the topic of the next section. 



III. DAMPING 

Thus far we have not considered dissipative or damp- 
ing effects tending to destroy the quantum coherence of 
the system; we have neglected the D term in Eq ||^ . The 
quantity Ptr has the interpretation of the degree of phase 
coherence between the two states, and D gives the rate 
of loss of this coherence 1^. In particular, we are inter- 
ested in the effect of D on the inversion process. With 
the loss of phase coherence between the two states, we 
expect the situation to become more "classical" and for 
the inversion to be inhibited. Indeed, in solving Eq[^ for 
large D one finds that the inversion is strongly blocked 
and that one arrives in the "Turing- Watched Pot-Zeno" 
regime 0. Eq 1^ with moving V(t) has been studied 
numerically Q. One finds that while for weak damping 
D has little effect on the inversion, values of D/Vtr ^0.1 
are enough to stall the inversion and values D/Vtr >> 1 
block it altogether. Given adiabatic conditions, the im- 
portant relation is that between the sweep speed and the 
decoherence time 1/D. This relation determines if the 
system has time to "decohere" during the sweep. For 
D/Vtr =0.1, for example, one finds that a factor of ten 
increase in the sweep speed will change the stalled in- 
version into a successful one. Such behavior is of great 
interest to us here, since if the experimental situation is 
favorable, we can pass from the quantum regime (small 
D) to the classical regime (large D) by varying the sweep 
speed. In the former case the system ends up in the other 
potential well while in the latter case, where the phase 
coherence is destroyed, the system behaves classically as 
Ptr ^ and remains in the original state. 

Concerning the absolute magnitude of D, a number 
of subtle issues are involved which cannot be treated 
here. However, for orientation we can take the result 
of calculations based on the Caldeira-Leggett approach, 
where we may identify the "decay rate" for weak dis- 
sipation with D. Weiss and Grabert [|l^ give T/Re^ 
for this parameter (their F). With R = 5 Mfi, this 
gives D = 0.08 mk = 9.6 MHZ at T = 100 mk, and 
D = 0.008 mk = 960 kHZ for T = 10 rriK. (Units: 
1 K ^ 8 meV = 120 GHZ, l/e^ = 4 k^). If these es- 
timates are correct, we are in an interesting range since 
it is possible to suppose sweeps both slow and fast on 
these scales and thus the direct measurement of D and 
its temperature dependence. 



IV. TIME SCALES 

Since we envision working at time scales shorter than 
have been common in this field, it may be useful to give 
a qualitative discussion of the various time scales which 
arise. The highest frequency involved is the ordinary har- 
monic frequency loq in one of the potential wells, giving 
the approximate level spacing in the well. For typical 
conditions this spacing may be on the order of several 
Kelvin (K). Since we envision experiments in the Kelvin 
to milliKelvin range, this implies that for the equilib- 
rium system at the start of the sweep only the lowest 
level is populated. The next parameter is the tunnel- 
ing frequency through the barrier. This is a sensitive 
function of the SQUID parameters, with sample values 
[/3l = 1.8, C = .01 pF,L = l Nh, i? = 5 MQ] we obtain 
« 20 GHZ, and ujtunnei/t^o ~ 4 ■ 10"^ 

The next two parameters concern the extent and speed 
of the sweep. The initial asymmetry e (Fig. 1) should be 
small compared to ujq in order to retain the approximate 
two-state character of the system, but large compared 
to ujtunnei to avoid initial mixing of the two states, say 
e/wo ''-^ 10~^ — 10~^. The speed of the sweep, ujsweep will 
be the easiest experimental parameter to control, and 
the behavior of the results as cusweep is varied will be an 
important check on the theory. It must not be so fast as 
to lose adiabaticity, but not so slow as to allow relaxation 
processes to mask the results. We may suppose it to be in 
the range lO-^t^o- 10"*c^o = 10-100 MHZ. With these 
numbers, it is possible to satisfy the adiabatic condition 
Eqp. 

Finally there are dissipative parameters D, and tOreiaxi 
the relaxation rate for transitions among the SQUID lev- 
els. From the above estimate of 1/D in the 1-10 MHZ 
range we have, as mentioned, the interesting possiblity 
of being able to choose sweep speeds either slow or fast 
with respect to the decoherence time while still retaining 
the adiabatic condition. The resulting switches between 
classical and quantum behavior would provide persuasive 
evidence for the correctness of our general picture, and 
allow the measurement of this parameter and its tem- 
perature dependence. Finally, the relaxation parameter 
^reiax givcs the rate of conventional kinetic relaxation 
processes induced by the environment, say thermal jump- 
ing over the barrier or direct transitions with the emis- 
sion of some energy. Estimates based on the formulas of 
Ovichinnikov and Schmid pl| ] lead to ujreiax ^ 100 kHZ 
which, as desired, would be much slower than the sweep 
time. This large difference between the tunneling time 
and the relaxation time is due to the many orders of 
magnitude suppression associated with the emission of 
energy in the relaxation process, particularly for under- 
damped SQUIDS. 

In conclusion, there are five important (inverse) 
timescales; ujq, uJtunnei, e, (^sweep, D, and ujreiax- Nat- 
urally a wide range of behavior is available by changing 
the parameters and some, like the sweep time or the tem- 
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perature, can be adjusted online. Probably the most im- 
portant condition for feasibility is the smallness of ujreiax '■ 
relaxation processes should not be significant during the 
sweep and the subsequent readout. 

Concerning the readout itself, we have examined a 
scheme involving a switchable flux linkage to a DC 
SQUID, which would be sufficiently fast for the above 
estimates . Here one profits from the fact that obser- 
vation of the system is only necessary after the procedure 
is completed. 



V. QBITS AND THE QUANTUM COMPUTER 

The two-state system under discussion here suggests 
itself as a physical embodiment of the "quantum com- 
puter" . The "qbit" itself is naturally represented by L 
and R playing the role of and 1. A linear combina- 
tion of L and R may be created by adiabatically rotating 
P from some starting position. Adiabatic inversion is 
evidently an embodiment of NOT since it will turn one 
linear combination into another one with the weights of 
L and R interchanged. As for CNOT, the other basic 
operation, a NOT is performed or not performed on a 
"target bit" according to the state of a second, "control 
bit" , which itself does not change its state. One straight- 
forward realization of this would be to perform the NOT 
operation as just described in the presence of an addi- 
tional linking flux supplied by a second SQUID nearby. 
The magnitude and direction of this linking flux would 
be so arranged that the inversion of the first SQUID is 
or is not successful depending on the state of the second 
SQUID. This and many other interesting combinations 
of junctions and flux linkages may be contemplated and 
are under study p^ . 

SQUID systems like these would seem to be partic- 
ularly well suited for the embodiment of the quantum 
computer, where we wish to generate a series of uni- 
tary transformations for the various steps of computa- 
tion. This may be done by creating a "moving land- 
scape" of potential maxima and minima, as in our sim- 
plest one-dimensional example of the adiabatic inversion. 
This imaginary landscape can be produced and manipu- 
lated by controlling various external parameters (as with 
our sweeping flux), performing the various operations in 
one physical device. Naturally, practicality will depend 
very much on the relation between the speed of these op- 
erations and the decoherence/ relaxation times which we 
hope to determine by the present methods. 

Although our interest here has been the SQUID, it will 
be evident that the principle of determining decoherence 
times through the inhibition of adiabatic inversion could 
be applied to many other types of systems as well. 

Preliminary reports on this work were presented by 
one of us (P.S.) at ICSSUR VI (NapoH, Italy, 24-29 May 
1999), at NATO ASI ( Ankara, 13-25 June 1999), and at 
MICRO AND NANO CRYOCENICS (Jyvaskyla, Fin- 



land, August 1-3 1999). 
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